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ABSTRACT 

Early reionization changes the pattern of anisotropies expected in the cosmic mi- 
crowave background. To explore these changes, we derive from first principles the equations 
governing anisotropies, focusing on the interactions of photons with electrons. Vishniac 
(1987) claimed that second order terms can be large in a re-ionized Universe, so we derive 
equations correct to second order in the perturbations. There are many more second order 
terms than were considered by Vishniac. To understand the basic physics involved, we 
present a simple analytic approximation to the first order equation. Then turning to the 
second order equation, we show that the Vishniac term is indeed the only important one. 
We also present numerical results for a variety of ionization histories [in a standard cold 
dark matter Universe] and show quantitatively how the signal in several experiments de- 
pends on the ionization history. The most pronounced indication of a re-ionized Universe 
would be seen in very small scale experiments; the expected signal in the Owens Valley 
experiment is smaller by a factor of order ten if the last scattering surface is at a redshift 
z ^ 100 as it would be if the Universe were re-ionized very early. On slightly larger scales, 
the expected signal in a re-ionized Universe is smaller than it would be with standard 
recombination, but only by a factor of two or so. The signal is even smaller in these 
experiments in the intermediate case where some photons last scattered at the standard 
recombination epoch. 
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1. Introduction 



Measurements of the cosmic microwave background (CMB) are very powerful probes 
of theories of structure formation. Even the so-called "medium-angle" CMB experiments 
arc looking at scales larger than those in the largest redshift surveys. This means that 
CMB observations are some of the purest observations in the field: more than any other 
type of experiment they sample the primordial spectrum of fluctuations. We can therefore 
hope to use the CMB experiments to answer many questions about primordial fluctuations: 
Are they Gaussian? What is the spectral index? Are they adiabatic or isocurvature? Do 
gravity waves play any role? 

These are profound questions and the possibility that we might answer them is one 
of the great promises of observational cosmology. However, there are several flics in the 
ointment. Perhaps most disturbing is foreground: our galaxy has lots of dust, synchrotron 
and free-free emission. All of these can get in the way of detecting the purely cosmic 
radiation. Another potential barrier to inferring the primordial spectrum from CMB mea- 
surements - and the one we wish to focus on here - is the possibility of reionization. In the 
standard cosmology the electrons and protons "recombinc" at redshift z ~ 1100, thereby 
cutting off contact between photons and matter. However, it is possible that hydrogen was 
"re-ionized" at a later epoch. The subsequent contact between photons and free electrons 
changes the pattern of anisotropics in the CMB. Thus we can use the CMB to infer the 
ionization history of our Universe. As we have said, knowing this history is crucial for 
the sake of determining the primordial spectrum, but it is also interesting in of itself. For 
example, if we determine that the Universe was ionized at a redshift z ~ 100, we will have 
learned something very useful about structure at that epoch. 

In this paper, we would like to quantify these statements with particular attention to 
experiments. The simplest way to discuss anisotropics in the CMB is in /— space (Bond 
et al. 1991). That is, we expand 5T{d, (j))/T = '^i^ cbimXimi^i (f>) define Ci = (|a/^p), 
where the angular brackets denote ensemble averages. Each theory has its own set of 
predicted Q's and therefore predicts that a given experiment will observe a variance 



where Wi is the window function appropriate to that experiment. 

FIG. 1. Window functions for several experiments currently searching for anisotropics in the cosmic 
microwave background. The curve labeled SP91 denotes the filter for the ACME telescope for the nine- 
point scan of 1991 (Gaier et al. 1992) in the South Pole; TENERIFE the one used by Davies et al. (1987) 
at Tenerife; MAX the filter used by Meinhold et al. (1993) in the region of /X-Pegasus; OVRO the filter 
used at Owens Valley (Readhead et al. 1989). For reference the COBE (Smoot et al. 1992) filter is also 
plotted. 

Figure 1 shows the window functions of a variety of experiments. For large I, the 
expected variance in a given experiment is seen from Eq. (1-1) to be roughly 1/2tt 
J (iln(/) {PCi)Wi. So the quantity ipCi) is convolved with the window function to give 
the expected variance. Figure 2 plots the f'Ci predicted by cold dark matter (CDM) with 
standard ionization history. The peak at I = 200 is probed by several of the experiments 
shown in Figure 1. 
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FIG. 2. The moments predicted by CDM with standard recombination. The quantity plotted l{l -\- 
X)Ci /QC2 is equal to one at low I where only the Sachs-Wolfe effect is significant. In the post-COBE era, 
we know that l{l + 1)Q/6C2 = + l)Q/(2.4 X 10-^)2. 

We want to know how the Cis change if the Universe is re-ionized. Is there still a peak 
at / = 200? Does the amplitude change? Does the peak move, perhaps to lower or higher 
11 Or both? To answer these questions we derive the fundamental Boltzmann equation 
governing the interaction between photons and electrons, the interaction that, as we will 
see, is responsible for the peak in Fig. 2. This derivation is presented in sections in what wc 
think is a very systematic way. One of the advantages of this systematic treatment is it will 
enable us to pick up not only the linear terms [which of course have already been derived 
many times (Peebles & Yu 1970; Wilson & Silk 1981)] but also the second order terms. 
These are of particular interest because Vishniac (1987) has analyzed one such second 
order term and concluded that in a reionized Universe, anisotropies are generated at small 
scales [high I]. By introducing this systematic treatment of the Boltzmann equation, we 
will be able to see if there are other second order terms that are as large as those considered 
by Vishniac, and if cancellations occur. 

First though, in section 4, we solve the first order equations. The focus is on two possi- 
ble ionization histories: the standard scenario wherein electrons and protons recombine at 
2; ~ 1100 with no further ionization and the opposite extreme wherein electrons remain free 
throughout the whole history of the Universe. In addition to presenting the results of a full 
numerical treatment, we also solve the equation [or an approximation to it] analytically to 
gain insight into the location of the peak in I— space. We find that this peak shifts to lower 
I if the Universe is re-ionized early enough. If the Universe never recombines, then there 
is a peak at I ~ 50. The signal in an experiment with a filter centered around this value of 
I would see a larger signal in a completely re-ionized Universe than in one with standard 
recombination. We also show why I'^Ci falls off at high /. This will lead into a discussion of 
the second order terms expected to be important. Sections 5 and 6 use this information to 
complete the derivation of the second order Boltzmann equation and present the solution. 
Our techniques are similar to those of Efstathiou (1988). In particular we will recapture 
Vishniac's result and show that the term he identified is indeed the dominant one. Thus 
we support Vishniac's conclusions about second order effects. Finally in section 7, we get 
more specific about experiments, presenting the predicted ((AT/T)g^p^.) in a variety of 
experiments for a variety of ionization histories. 

One final introductory comment: The signal from the cosmic background depends on 
many parameters even in a model as simple as CDM. For example, the spectral index n 
need not be equal to one, and even a slight deviation has dramatic implications for the small 
and medium scale measurements. The signal in a given experiment also depends (Bond 
et al. 1991; Dodelson & Jubas 1993) on Q^, the fraction of critical density in baryons 
today, and h which parametrizes the Hubble constant today {Hq — 100/ikmsec~^Mpc~^). 
And of course, models other than CDM give different predictions and have different sets of 
parameters to fiddle with. Here we are focusing only on the impact of different ionization 
histories. It therefore makes sense to fix all these other parameters and play with the one 
variable of interest: the ionization history. Accordingly, we will focus only on a cold dark 
matter dominated Universe with Harrison- Zel'dovich spectrum {n = 1) and set h = 1/2 
and Qb — 0.05. Several other groups (Bond & Efstathiou 1987; Efstathiou 1988; Hu, Scott, 
& Silk 1993) have recently considered the effects of re-ionization for other cosmologies, in 
particular the minimal isocurvature model proposed by Peebles (1987). 
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2. Compton Collision Term: A General Derivation 

The photon spectrum is governed by the Boltzmann equation: 

^/(x,p,t) = C(x,p,t). (2.1) 

Here / is the photon occupation number, a function of momentum p, position x, and 
time t. The coUision term, C, also depends on these variables. Theoretically, it includes 
contributions from all scattering processes, although in practice only Compton scattering 
off free electrons need be considered. 

In the absence of collisions [C = 0], Eq. (2.1) says simply that photons travel freely 
along geodesies. For example, in a Robertson- Walker background, the left hand side of 
Eq. (2.1) becomes 

where H is the Hubble rate and p = |p|. If we are interested in large scale anisotropics, we 
must do better: we must account for the perturbed metric when expanding d/dt in terms 
of partial derivatives. Such an account leads to the Sachs-Wolfe effect (Sachs & Wolfe 
1967). [Recently, Martinez-Gonzalez, etal. (1992) have taken this program a step further 
and considered second order effects in the perturbations to the metric] 

In this paper we will focus on the right hand side of Eq. (2.1): the collision term. This 
term governs small scale anisotropics and spectral distortions. In the limit of completely 
elastic collisions, the right hand side vanishes. Typically, in the regime of interest, very 
little energy is imparted from electrons to photons in a collision, so this limit is a good 
approximation. To get non-zero effects, therefore, we simply need to expand the right 
hand side systematically in powers of the energy transfer. 

The starting point then is the collision term corresponding to the process 

e(q)7(p) ^ e(q')7(pO- 

To calculate this, we find the matrix element M for the process, square it, weight it by 
the occupation numbers of the particles and integrate over all other momenta, q, q',p'. 
Therefore, 

^^^ = 11 (2.)l2l(q) (2.g(qO (2.g(pO P-''^^'^ ^ " - 
S(q')/(P') (l + /(P)) - S(q)/(p)(l+/(p')) 



X 



(2.3) 



where E{q) = \/ -\- w? and the delta function enforces energy- momentum conservation. 
The last line in Eq. (2.3) contains the distribution functions, / that of the photons and g 
the electrons. We have dropped the Pauli suppression factors 1 — g-, since in all realistic 
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cosmological scenarios, g is very small. We don't know / [that is what we are trying to 
solve for] but we do know g: Due to the fast rate of Coulomb collisions, the electrons are 
kept in thermal equilibrium, so 



^(q) = rie 



27r 



3/2 



exp 



-(q- mev) 



(2.4) 



where v is the velocity of the electrons and the normalization comes from requiring that 
J d^q g/{27r)^ = rig, the electron density. Before beginning the expansion, we can trivially 
do the q' integration in Eq. (2.3) by using the three dimensional delta function. This 
leaves 

C(p) = dp'p'— / — -— J — ^ -6{p + E(ci)-p' -E(p + ci-p')] 



X 



^(p + q-pO/(pO 1 + /(P) 



^(q)/(p) l + /(p) 



(2.5) 



As mentioned above, the energy transfer, -E'(q) — E{q + p — p'), is small compared 
with the typical thermal energies which are of order T. In fact, the energy difference is of 
order £'(q) — £'(q+p — p') ^ (p' — p)-q/TOe = 0{Tq/me). Thus our expansion parameter, 
the energy difference over the temperature, is actually q/rrif.. The electron momentum has 
two sources: the bulk velocity {q = niev) and the thermal motion (q ~ ^/m^). Thus, 
an expansion in q/rrie is necessarily an expansion in v and ^/Tjrnl. At the end of our 
expansion we will have a first order source term linear in v and a host of second order 
terms quadratic in v and \jT/me- 

The strategy now is to expand everything - energies, squared matrix element, delta 
function and distribution functions - using the energy transfer as an expansion parameter. 
The Boltzmann distribution expands to 



5r(p + q-p') =^(q)h- 



(p - p') ■ (q - mv) (p-p 



l\2 



1 

+ 2 



(p - p') ■ (q - ^v) 

TUeTe, 



+ . . . 



(2.6) 



The second term in brackets is first order in the perturbative quantities while the last two 
are second order. Meanwhile the delta function expands to 



5{p + E{ci)-p' -E{p + q^-p')] =5{p-p') + 



(p-p') ■(id5{p-p') 



rue 



dp' 



{p - p'f d5{p - p') 1 
2me dpf 2 



(p - p') q 



rup 



2a2 



d'^5{p-p') 



dp 



n 



+ 



(2.7) 
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The derivatives of the Dirac delta functions here wiU uhiniatcly be handled by integration 
by parts. [Bernstein (1988) introduced this approach to derive the Kompaneets equation 
which, as will see in section 5, is a special case of the general second order equation.] 
Finally, we expand the photon distribution function. 



(2.8) 



where (p) is the zero order photon distribution [typically Planckian] which of course 
depends only on the magnitude of p. We can use these three expansions to rewrite Eq. 
(2.5) as 



^^^^ - 8^ y ^ 4^ y (27r)3 E(q)£;(p + q - p') ^ 



) + 



(P-P') -f^d^P-p') 



rUe 



dp' 



+ 



(P 



- Y>'f dd{p - p') 1 
2me dp' 2 



(P - P') ■ q 



n2 o2 



d S(p — p') 



dp 



n 



X 



+ 



j^^\jj)-r\p) 



(0), 



+ 



/^^Hp')-/^'np)-/^'np')(i + /^'np)) 



(p-pO • (q-mev)" 



TTlp l p ^ ' 



+ f^'^ iP'Ki + f^'^ iP)) ^ 1 ^ (P - pO • (q - -ev) ^ 2 



2mpTp. 



rUeTe. 



) 



(2.9) 



This equation looks like a mess, and we have not even expanded all the quantities on the 
first line yet. It turns out though that the hard part is over. We now recognize that 
the zero order term in Eq. (2.9), i.e. the one we get when multiplying together the first 
terms in each of the curly brackets, vanishes. That is, 5{p — p'){f^'^\p) — f^^\p')) is zero. 
Therefore only terms of first order remain after multiplying all the terms in the two curly 
brackets. This means that [since we are interested only in terms up to second order] we 
only have to keep first order terms when we expand the matrix element and the energies 
on the first line. In fact to first order the energies can be simply replaced by rrie, and the 
matrix element squared is 



|Mp = Gnarml ((1 + cos26') - 2cos6'(l - cos6')q • (p + p')/me + . . .) 



(2.10) 



where cos^ = P • P' and ctt is the Thomson cross section. To simplify things further, 

we can explicitly do the d^q integral by using {g) = ((/q) = neTOgV; and {gqiQj) = 
SijUerrieTf, + riemlviVj, where (A) = J d^q A/{2it)^. Then we find 



C(p) = / dpp 



4p 



cW(p,pO+c5J)(p,pO + 4'^(p,p') 

+ cS(P,P') + 4inP,P') 



(2.11) 
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where the first order integrand is 



c(^)(p, p') = (1 + cos^^) Sip-p')if^'\p') - /(^)(p)) 



+ {f^'\p')-f^'\pmp-p')--^^^^y, (2.12) 

and we have separated the second order terms into four parts. The first set of these 
contribute to what is referred to as the Kompaneets equation (Kompaneets 1957) describing 
spectral distortions to the CMB (Bernstein & Dodelson 1990): 



4\p,p') = {l + cos'9) 



/\2 



2M (P - P') 



2mp 



dp 



n 



[f('\p') + f('\p) + 2f('\p')f('\p)) 



dS{p — p') 
dp' 



2{p - p')cose{l - cos^d) 



rUf 



6{p-p')f^'\p')il + f^'\p)) 



- {f^'\p') - f'\p) 



dS{p — p') 
dp' 



There is also the simple damping term 

4'^(p,p') = (1 + cos^^) 5{p - p'){f^^\p') - f^'\p)); 
a set of terms coupling the photon perturbation to the velocity 



(2.13) 



(2.14) 



.(2) 

"Ad 



(p,p') = (i+cos2^) (/«(po-/^'np; 



(l + cos^^) (p-p')-v-^'^(^-^') 



dp' 



2cos6'(l - cose)5{p - p'){p + p') • V 



(2.15) 



and finally a set of source terms quadratic in the velocity: 
cl?(p,p') = (/^°^(p')-/^°np)) (P-P')-V 



(1+cos^^) ip-pr-d^s{p-p') 



dp 



/2 



2cos6'(l-cos6')(p + p')-v 



d5{p-p') 
dp' 



(2.16) 
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We'll see in the next section that the first order terms, those in Eq. (2.12), reduce to the 
standard first order equation once the and Q' integrals are done. There are a lot of 
second order terms. As we have mentioned, one subset of these terms has already been 
studied: those leading to the Kompaneets equation (Kompaneets 1957). We know that 
Vishniac analyzed a second order term. Which one is it here? It turns out that he looked 
at none of the second order terms in c*^^-* [or at least he didn't write a paper about any of 
them]. Rather, he got yet another second order term by expanding the electron density as 
TT-e = ?ie(l + ^^e) and then multiplying de by the first order terms in Eq. (2.12). 

3. First Order Equation 

In this section, we derive the well-known first order equation coupling photons and 
electrons. To do this, we need focus only on the terms in Eq. (2.12). For the angular 
integrals, we choose the polar axis to lie along the direction of the electron velocity, so that 
azimuthal symmetry is maintained. Then //' is the polar angle defined by //' = v • p'; we 
also define = v • p. Thus we have 

^^'np) = ^ /^^>V^C«(P,P') 

jyp'p'\s{p-p') 1'^^ (/W(pO-/W(p))j^''^^(i + (p.pO^) 



Ap 



+ «(/'»>(?') - /<"' W)^^ £ f (P^ - pV) f (1 + (P ■ P')=) 



(3.1) 



where v = |v|. The dot product p-p' is a function of both ji' and (/>'. But /'^^-'(p') depends 
only on p' and ^' , not on (j)' because of the azimuthal symmetry. [This is perhaps too 
strong a statement. For many metrics considered by cosmologists, there is an azimuthal 
symmetry so /*^^-* depends only on the polar angle. There are cosmologies though wherein 
this symmetry is not maintained, so f'^^^ could well depend on the azimuthal angle. This 
happens for example when tensor perturbations (Crittenden et al. 1993) or cosmic strings 
are present. We will restrict our analysis to cases where the symmetry exists.] To do the 
(j)' integral we first rewrite the integrand in terms of Legendre polynomials 

+(P-P')' = ^(l + ^^'2(P-P')) (3.12) 
By the addition theorem of spherical harmonics, 

J^2(p-pO= E W-^^2"^(P-v)P2'"(P'-v)e^™^<^'-'^^ (3.13) 



m=-2 
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Thus, 



/ ^P2ip ■ pO = P2(P • V)P2(P' • V) ^ P2(/X)P2(/X') (3.14) 

To do the //' integrals, it is useful to define the moments of the distribution function: 



fiip)= t ^Piil^)fip,f^). 



(3.15) 



-1 



Now we can use the orthonormality of the Legendre polynomials to write 

' ^(l + i^2(/i)P2(/x')) [f^'\p')-f^'\p)]-f!,'\p') + lA'\p')P2{^i)-f^'Hp,f^) 

(3.16) 

The notation at this stage is a bit confusing, so let's restate it: the superscript refers to 
the order of perturbation theory. Here we are considering the first order correction. The 
subscript refers to the moment of the distribution function; note that the last term in Eq. 
(3.16) has not been integrated over, i.e. still depends on so it has no subscript. Eq. 
(3.1) now becomes 

c^'\p)='^ dp'p' sip-p') (/r(pO + ^/i'^(p')^2(/.)-/«(p,/x)) 



+ vp^{f'\p')-f^'\p))^^^^ 



(3.17) 



The remaining p' integral can be done, in the first case trivially and in the second by 
integrating by parts. Thus our final expression for the collision term is 



C^^\p,IJ,) = UeCTT 



1 



_ f(l) _ 



a/(o) 
dp 



p- 



(3.18) 



A(/x) 



(3.19) 



It appears as if this collision term is momentum dependent. To get rid of this illusion, it 
is useful to define 

-p dfioy 

4 dp 

Then combining the left hand side of the Boltzmann equation from Eq. (2.2) and the right 
hand from Eq. (3.18), the full first order equation is 

'''''^ + ^ A2P2(//) + ^fiv - (3.20) 

[The expansion term in Eq. (2.2) simply forces the zero order distribution function to 
depend on the comoving momentum pa where a is the scale factor. One can then show 
that the expansion term drops out of the first order equation for A.] The subscripts here 
again refer to the moments of A(/x) defined just as in Eq. (3.15). Since to first order A 
is independent of photon energy, p, A as defined by Eq. (3.19) is equal to the brightness 
defined as 

^"^■'f"' (3.21) 



A = 
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4. First Order Solutions 

Now that we have the equations describing the interactions of photons with electrons, 
we can solve them to determine the predicted anisotropics in the cosmic microwave back- 
ground. Figure 2 shows the results of numerically integrating the full set of linear equations 
starting with these initial conditions, assuming standard recombination at 2 ^ 1100. We'd 
like to do two things in this section: First, we would like to understand the peak at / ~ 200. 
If we understand why it occurs at / ~ 200 when the Universe follows a standard ionization 
history, then we will be able to understand how this peak shifts when we consider differ- 
ent ionization histories. Second, we would like to understand the damping that occurs at 
I ~ 1000. We will see that this damping is due to the "finite thickness of the last scattering 
surface." For a re-ionized Universe, this thickness is much larger, and therefore, damping 
is apparent even on very large scales [much smaller /] . It turns out that understanding the 
physical reason for this damping will give us a clue as to which second order terms are 
likely to be significant. [Due to the crudeness of our approximations, we will not be able 
to account for the oscillations that occur as Ci is damped; these are due to acoustic waves 
at the time of recombination.] 

To understand these features, we will solve a simplified version of Eq. (3.20). First 
note that A depends on position; to account for this, we can Fourier transform it: 

A(x, p,t) = V J (fk A{k, r)e^^p'^^-"'-^, (4.1) 

where /cphysicai — k/a; k being the comoving wavenumber; a the cosmic scale factor; V the 
volume, a factor which drops out of all physical results. Note that because of the azimuthal 
symmetry, A depends only on the magnitude of k and the dot product k • p = v • p = /x. 
Note that the first equality holds here since the velocity is irrotational, meaning that the 
Fourier transform of the velocity is parallel to k. We also define the conformal time 




The conformal time today is 

2 

To = — = 1.2 X 10^^ sec (4.3) 
-no 

since we're assuming a flat, matter, dominated Universe with h = 0.5. The first order 
equation is now 

(A/iv - a) (4.4) 

where the dot denotes derivative with respect to r. The Ao and A2 terms on the right 
hand side of Eq. (3.20) have been neglected. That's because we are only interested in 
the kick that photons get from the moving electrons [the A/iv term] , the so-called Doppler 
effect. The Aq, A2 terms give rise to what's sometimes called the intrinsic anisotropy. Also 
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we still have not written down the terms that represent the perturbation to the metric. 
These lead to the Sachs- Wolfe effect. So all we should get out of the simplified Eq. (4.4) 
is the Dopplcr effect. 

Eq. (4.4) is a first order differential equation, whose solution is 

A(/c, II, t) =A(/c, jjL, Tinitiai) exp I i/c//(Tinitiai - t) - dr' ne{T')aTa{T') \ 

^ "'Tinitiai J 

+ 4// /" dT'ne{T')(TTa{T')v{k,T')e'^^^'''-^^ 

•'Tinitiai 

X exp |- dT"ne{T")(7Ta{T")^ . (4.5) 

The first term on the right represents the anisotropies that were initially present and have 
persisted until some late time r. These are damped out by scattering with electrons, the 
/ neOTddr' term in the exponential. If we choose an early enough rinitiab then this integral 
is always large, and so the initial anisotropies are not important today. We can write Eq. 
(4.5) in a more compact form by setting Tinitiai to zero and by defining the visibility function 

g{T, t') = ne{T')aTa{T') exp |- dTW{T")<TTa{r")^ . (4.6) 

Then we have simply 

A(A;,//,t) = 4// / dT'v{k,T')g{T,T')e'''^^^'-^\ (4.7) 
Jo 

The visibility function defined in Eq. (4.6) has an interesting physical interpretation. 
To see this, first note that dT'g{oo, t') = 1, so ^ is normalized like a probability density. 
In fact, that's what it is: the probability per unit conformal time that a photon at time r 
was last scattered at time r' . The visibility function depends on the free electron density 
Ue] thus it is particularly sensitive to the ionization history. Figure 3 shows the visibility 
function for two ionization histories: (i) standard recombination at ^ ~ 1100 and (ii) the 
case where the electrons remain ionized throughout the whole history of the Universe. In 
the first case recombination happens rapidly, so the "surface of last scattering" is centered 
at ^ = 1100 but with a very small width. In the second case, the surface of last scattering 
is centered at 2; = 100 but is quite wide. To simplify Eq. (4.7) further, we can approximate 
5r as a Gaussian: 



1 r {r'-TR 



i2 



giro, t') ~ exp <^ ^ } . (4.8) 



2 

'R 



For standard recombination, tr ~ .02to while the width of the last scattering surface, 
Str ~ -Itr. For no recombination, tr ~ .OStq while Str ~ .06to. 
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FIG. 3. The visibility functions for standard recombination [solid line] and no recombination [dashed line]. 



Our goal is to find out how A depends on the ionization parameters, tr and Str. 
With the approximation in Eq. (4.8) we can go further and solve explicitly for A in Eq. 
(4.7). We now have 



A(/c, yu, To) = 4/x / dT'v{k, r')e^^'^^^'-^°) exp <^ ■ 



jr' rg) 
8rl 



(4.9) 



In a matter dominated Universe, v grows as r, so v{t) = v{tq)t/tq. We can also set the 
upper limit of the integral to infinity since the exponential is negligible for t > tq. Then, 



^/-KTQdTR I dikfj.) Jq 



5rl 



(4.10) 



where we have written t' as {l/i)d/ dikfj). The integral in Eq. (4.10) contains a lot of the 
physics we are after. 

FIG. 4. Real part of the integrand of Eq. (4.10). 

Figure 4 shows the integrand for two limiting cases: (i) kjidrn » 1 and (ii) kiiSr^ « 
1. In the first case [large fc], the scale of the perturbation is much smaller than Str, the 
width of the surface of last scattering. A photon travelling through the last scattering 
surface travels through many regions where v is positive but an almost equal number of 
regions where v is negative. Thus the total contribution to A is small on scales smaller 
than the thickness of the last scattering surface. The second case comes about either 
because k is small [i.e. the scale of the perturbation is large] or ~ [p • v oc p • V5 ~ 
0, the photon is travelling perpendicular to the direction in which the perturbation is 
changing]. Perturbations on large scales do make a large contribution to A since there is 
no cancellation through the last scattering surface. Similarly, there is no cancellation [in 
the integral] if the photon is travelling perpendicular to the gradient of the perturbation. 
[Note though that Eq. (4.10) has a factor of /x in front which does lead to a cancellation 
for such photons.] In any event, the integral under discussion can be written down in 
terms of Error functions, but for our purposes we will make another approximation which 
will simplify things further. The contribution of the integrand at the lower limit r' = is 
suppressed by a factor of e~^'^^/^'^^^ ; for most realistic ionization histories this will be a 
pretty small number, so we can extend the lower limit all the way to — oo with little loss 
of accuracy. The remaining integral is then the Fourier transform of a Gaussian, which is 
itself a Gaussian, or 



AD°PPi«'^(A;,//,To) = 4//i}(A:,To) 



To 



1 + i 



kiidr^ 
2tr 



giA;/i(rH-To)g-(fep(5rH/2)^ 



(4.11) 
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where the superscript indicates that this is the perturbation to the photons caused by 
the Doppler effect. Eq. (4.11) clearly illustrates the damping on small scales as we have 
discussed; the damping scale is of order k ~ 1/Stji. Kaiser (1984) pointed out that since 
the velocity is parallel to k, there is an extra factor of n in front of Eq. (4.11). Thus, 
for example, dp^/p^ = (1/2) J\dp A^^PP^'^'ik, p,To) oc l/{k6TR)^. This suppression of 
Sp^/ Pj would be averted if the source term had a component perpendicular to k. The other 
feature of Eq. (4.11) which is of interest is the factor of tr/tq in front. This tells us that the 
later the photons are in contact with the electrons, the greater is the Doppler kick they get. 
But this is what we expect: velocities grow with r, so photons scattering off electrons at 
later times are seeing larger velocities. This simple fact will have surprising consequences, 
namely, on large scales [before damping sets in] a reionized Universe produces a larger 
signal in the CMB than one which never reionizes! 

What is the contribution of the Doppler effect to the present day anisotropy as mea- 
sured by the Q's? Referring back to Figure 2, we see that for Z < 30 or so, the dominant 
contribution to C; comes from the Sachs-Wolfe effect, for which /(/ + 1)C; is constant. 
At larger /'s, the Doppler effect becomes relevant. The reason for deriving an analytic 
expression for A^^pp'*^"" in terms of 6tr and tr is to see when the Doppler effect becomes 
important as we vary the ionization history of the Universe. To calculate the Q's, we 
expand A^°PP'®''(fc, p, tq) in a series of Legendre polynomials, defining 



(4.12) 



Then, C/ is given by: 



^Doppler ^ / { A^'°PP'"'^(A;,To) 



V_ 
Stt 



dk k^ { 



dpPi{p)4pv{k,To) 



Tr . kpdr'^ 



To 2to 



R 



^ ^ikiJ,{TR-To)^-{kiJ,STii/2)^ 



(4.13) 



Here we are using the normalization of Efstathiou, Bond, & White (1992). [A good way 
to check normalization is to insure that C2 = (S/tvtq) dkj2{k)P{k)/k'^, where P{k) 
the power spectrum is normalized to COBE, so on large scales is proportional to C2A;.] 
As a first approximation to (jj^'^PP^'^^ ^ we'll assume that only k < 2/6tr contribute to the 
k— integral due to damping, and that for these values of k we can set the exponential 
^-(kfiSTn/2) Qj^g Then the integral over p is simply the derivative of a spherical Bessel 
function: 

2 7-1 id{k{TR-To)) {-ly 



(4.14) 
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and similarly the term with fi^ becomes the second derivative of a spherical Bessel function. 
Also the continuity equation tells us that v{k,To) = —i'k5e{k,To)/k^ = —2i'k5e{k,To)/k'^To, 
where Sg is the fractional change in density of the electrons, assumed equal to that of the 
rest of the matter. Thus, the ensemble average is: 



< \d{k,To)\^ >-- 



(kro) 



4 ,f,, m2 4P(/e) 

< \5e{KTo)r >- — 



v{kToy 



(4.15) 



Here P{k) is the COBE normalized power spectrum;for CDM (Efstathiou et al. 1992), 
PcDM{k) = ^C2T^kT^{k) where C2 = ^T^iQuMs/Tof = 9.7 x IQ-^^ and T is the 
transfer function. 

We now have our final result for (^^^opp'^"" in terms of a simple one dimensional integral: 



^Dopple 



TT \ To 



H 2 



dz 



2tqTr dz'^ 



(4.17) 



FIG. 5. l[l + 1)C/ for standard recombination and no recombination. The results of a full numerical 
solution of the coupled Boltzmann equations [solid lines] are compared with the approximate solution given 
in Eq. (4.17) [dashed hnes]. 

Figure 5 shows this approximate result for c^°pp^^^ added to the Sachs- Wolfe result of 
^Sachs- Wolfe _ QC2/ 1(1 + 1) for Standard recombination and no recombination. Also shown 
are the results of numerically integrating the full set of Boltzmann equations [including 
perturbation to the metric, CDM and three species of massless neutrinos] for these two 
ionization histories. Apparently the approximate result of Eq. (4.16) gives a very good 
qualitative picture of the results we are after: (i) the Doppler peak is larger on large scales 
[small I] and (ii) the damping does take place at larger scales than in the case of standard 
recombination. Clearly, though, for an accurate quantitative analysis the equations must 
be solved numerically. One final quantitative point: Having emphasized the fact that Ci 
is larger in a no-recombination Universe at low I, we should point out that it's not that 
much larger. The difference never exceeds 30% for any I. Thus we'll see in Section 7 that 
the expected AT [which goes as the square root of Ci] for no-recombination is roughly the 
same on these scales as for standard recombination. All our work in this section has been 
on the linear terms in the Boltzmann equation. We now turn to second order terms. 
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5. Second Order Equation 

In this section we write down the second order equation which follows from our general 
Boltzmann treatment in Section 2. It turns out that there are numerous terms. We can 
immediately throw out one class of terms in Eq. (2.11) though. The terms coupling f^^^ 
and V, i.e. those in Eq. (2.15), are nominally second order. In practice, though, we have 
seen that /(^^ [or its equivalent A*^^)] is very small on small scales due to damping. So we 
we can neglect c^^ in Eq. (2.11). The remaining terms can be manipulated as were the 
first order terms in Section 3. The result for the second order collision term is 

+ ,2^^(^^ + l)+,V^(^/.^ + ^) 



J_d_ 
ml dp 



(5.1) 



The first term on the right hand side arises from writing the electron density as Ue = 
ne(l + Se) and then multiplying the first order collision terms (Eq. (3.18)) by 6e [note 
that we have dropped all terms with f^^^^s in them]. This is the so-called "Vishniac 
term" . The remaining terms arise from the second order expansion in the collision integral, 
Eq. (2.11). At second order, electron-photon scattering is not purely elastic. Energy is 
transferred between the electrons and photons, inducing spectral distortions. In the limit 
of no anisotropy only the terms on the last line, the Kompanects terms which give rise 
to spectral Sunyaev-Zel'dovich distortions, survive. Since the Kompaneets terms do not 
induce anisotropics, we shall neglect them. Finally, since we are interested in scales much 
smaller than the horizon, non-Newtonian gravitational effects can be ignored. 

It would be nice to define a A*^^^ just as we did A^^) in Eq. (3.19) so that the 
second order equation is momentum independent. This is impossible though, because the 
second order equation is momentum dependentl Specifically the v'^ terms in Eq. (5.1) 
are momentum dependent. Thus in general A^^) will depend on momentum. There are 
two ways to deal with this: Either we can go ahead and try to solve the full equation for 
A^^)(A;,p, ii), or we can integrate out the momentum dependence. Let us first take the easy 
way out, and later we will see if it is necessary to include the full momentum dependence. 
To do the integration we can define 

as in Eq. (3.21). Then integrating Eq. (5.1) leads to the second order equation: 



aA(2) 
dr 
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where we have dropped the Ag and terms on the right hand side, as these should 
be irrelevant on small scales. On physical grounds, we can argue that the only important 
second order term is the "Vishniac" term. The Vishniac term is proportional to SeV whereas 
the other terms are proportional to f^. By the continuity equation, though, v ~ Se/kr. 
Since kr is quite large (> 100) for the scales of interest, the Vishniac term should dominate. 

Are the other terms in Eq. (5.3) completely irrelevant? Not necessarily. Recall that 
we integrated out the frequency dependence of the velocity squared terms in Eq. (5.1). 
Even if the final value of < (AT/T)^ > induced by these terms was a factor of 100 smaller 
than the leading terms, we might be able to pick them out because of the frequency 
dependence. [A recent example of the process of separating out a signal from a much more 
powerful source with a different "spectral index," i.e. frequency dependence, can be found 
in Meinhold et al. (1993).] So we choose to keep these additional terms a little longer. 

6. Second Order Contribution to ^ 

In this section we will follow the treatment of Efstathiou (1988) in solving the second 
order equation. We can rewrite Eq. (5.3) as 

55,(k,T) + 5^,(k,T)- A(2)] (6.1) 

where the Vishniac source term is 

SsviK r) = 4p • ^ v(k', T)5e(k - k', r) 

k' 

= (:^)'^E^'^"e(k',ro)<je(k-k',ro) (6.2) 
and the source term quadratic in velocities is 

5„,(k,T) = 5^[7v(k',T)-v(k-k',T) + 15p-v(k',T)p-v(k-k',T)" 
k' 

-4(t/to)^ ^ 7k--(k-k0 + 15p ■ k-p • (k - kQ ~ ~ 

= 2^ \k-k'\^k'^ 5e(k,To)5e(k-k,To). 

(6.3) 

Here we have Fourier transformed from position space to momentum space. We have 
also used the facts that (i) the velocities are first order, and are related to the baryon 
density perturbations through the equation of continuity, and (ii) the time dependence of 
the perturbations is given by Se{T) = 'Je(T"o)(~)^ in a matter dominated Universe. 

The second order equation is of the same form as the first order equation, whose 
solution we wrote down in Eq. (4.7). By analogy, the solution to the second order equation 
is 

A(2)(A;,//,To) = e-"'^'^° [^^.(k, to)/i(/c//to) + 5,,(k, to)/2(A;//to)] (6.4) 



A + ik/iA^'^^ 
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where the time integrals are 

Jo ^T-q/ 

hikfiTo) = r dT'(-Yg{To,T')e"'^'^' (6.5) 

with g our old friend, the visibility function. We encountered an integral of this form 
[only the powers of t' differ] in section 4 when analyzing the first order equation. The 
main lesson we learned was that integrals of this type are strongly damped unless kfi is 
small. On small scales [large k] this means that the integrals are non-negligible only if 
~ 0. Thus the main contribution to A^^^ can be found by evaluating the Ssv{to) and 
SvviTo) at = 0. Before doing this, let's go back to Eq. (4.7). There the source term was 
proportional to fi, so the contribution of the linear source term was greatly suppressed. 
Vishniac's profound observation and the origin of the dominance of second order terms is 
that the second order terms Ssv and Syy do not vanish at = 0. 

For the very small angle experiments where second order effects are important, there 
is a simple formula for the C^'s: 

IV 

Ci = rf/i {\A{k = l/To, /x)|2) (6.6) 

where V is the volume which will drop out at the end of the calculation. To derive this, 
start with the small angle formula (Doroshkevich, Zel'dovich, & Sunyaev 1978): 

Cie,a) = V j ^ j di,{l-fi'){\A{k,fi,To)\')Jo{kRce{l-fi')) e^p[-{kR,af] 

(6.7) 

where 9 is the angle between the two observing direction; a is width of the Gaussian 
beam observing the temperature differences; and i?c is the comoving distance to the last 
scattering surface, Rc ~ tq in a flat universe with Zc ^ 1. Since (|A(A;, |U, ro)P) is highly 
peaked around /U = 0, we can set ^ = everywhere else in the integrand. Meanwhile the 
product of Jo and the exponential suppression is just the experimental window function 
referred to in Eq. (1.1). Changing variables to I = k/ro leads to 

<^(^'^) = ^r^/ ^^^'/ d^i {\A{k = l/ro,^i)\^)Wi,,,pt. (6.8) 

047r Tq Jq J_i 

Comparing this to the continuous form of Eq. (1.1) leads to Eq. (6.6). Let us now calculate 
the ensemble average (|A(A; = Z/to,/x)|^): 

(|A(A;,//,To)|2) =< \Ssy\^ > < SsvS:, > hl^+ < S^Sy, > Ilh 

+ < \Syy\^ > I/2I'. (6.9) 
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We argued above that Ssv/Syy < 10~^, so the last term here is suppressed by a factor 
of at least 10~^ relative to the first and we can safely drop it. The only terms with Sw 
that might be interesting are the cross terms, so we will keep these. Upon squaring S, we 
encounter double sums, say over k' and k"; thus we need the identity 

< 4*(k04*(k - k%{k'%{k - k") >= ^s^,^^„ + <5j,,+k„ J,] (6.10) 

where P is the power spectrum today and S, the Kronecker delta. We can use this expres- 
sion to expand Eq. (6.9) as 

f\A ft ■ k" 

mk,^^,To)f) = ^Yl ^(mk - k') + W',k] ^ 

k',k" 

X / P ■ 1/ P + 7k--(k-k0 + 15p • k-p • (k - kQ 1 

To go further we use the Kronecker deltas to get rid of the k" sum and change the k' sum 
into an integral via V J d^k' / {2t{)^ . Furthermore, due to the damping effect we 

can set p ■ k = /X = everywhere except in the time integrals (/i, 12)- Thus, 

mk, .,r„)P) = / P(k')P(k - k') p . k' (i, - ^) 

With no loss of generality we orient our coordinate system such that: 

k-p = 0; k-k' = kk'x; k' • p = k'^/{l - a;2)sin0 (6.13) 



Then 



{\A{k,^i,To)\') = y^^^^ dk' J ^dx d<j> P{k')P({k^ + k'^-2kk'xf/^) 



k'^/i 



1 1 



X /e V i — a;^ sm ( 



2- ^iri2 , 7{kk'x-k'^) - 15k'^{l-x'')sm^<P^ 



k'^ A;2 + /c'2 _ 2kk'x 
-_k^ 

(^ + - 2kk'x)TQ 



(6.14) 
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and 



We now see that the cross terms - those with /1/2 - vanish once the integral is done (Hu 
et al. 1993). Thus we conclude that the only second order term of any significance is the 
Vishniac term. Plugging this expression into Eq. (6.6) yields 

lJP\k^l/To) K{k = l/To) 

= 2^1 (^-^^^ 

where 

Kik)~' Hdk' f-\ /(^'H^^' + ^"''^^'^)'^') il-.^)ik^-2.kk') 

J=^ j'^df,\h{kf,To)\'. (6.17) 

For CDM, the integral over the power spectrum, can be approximated by K{1) ~ 
-0.1 + 1.1(//1000) for 1000 < / < 5000, with a maximum error of 5%. [The integral K is 
half of the integral Efstathiou (1988) calls 12- It's also much easier to compute numerically. 
To prove the identity, one uses the invariance of the integrand under k' -(-^ k — k'.] 

FIG. 6. C^'s induced by the second order Vishniac term for CDM with no recombination. 

The integral J contains all the information about the ionization history, but it also 
appears to depend on wavenumber k. Efstathiou (1988) made the intriguing observation 
that for most reasonable ionization histories, J is in fact independent of k. To see this, we 
note that the integrand is sharply peaked around |U = so there should be very little error 
introduced if we extend the limits of integration all the way to ±cxo. Then, 

J^:^ dikfiTo) / dr' (-) 9{To,T')e'''^^ / dr" (-) <7(to, r'^e-*^'^^ 
J-00 Jo ^^0^ Jo ^'To^ 

= rfT'(-) /(ro,r')ro (6.18) 

where the last equality follows immediately since the (kfiTo) integral yields a delta function 
in t' — t" . All the information about the ionization history is in the integral J. From 
Eq. (6.18), we see that the integrand is heavily towards late times by the r'^ factor. 
Thus ionization histories wherein the visibility function is peaked at late times - that 
is, scenarios in which the Universe is re-ionized - are most likely to produce appreciable 
secondary anisotropics. For standard recombination, J = 1.7 x 10~^ while J = 7.3 x 10~^ 
for no recombination. 

Figure 6 shows the Q's for CDM with no recombination. The shape of the curve is 
exactly the same for other ionization histories, only the amplitude, which is determined 
by J, drops. Nonetheless, since the integrand of J is heavily weighted towards late times, 
even relatively late re-ionization produces a comparable signal. The other feature of note 
in Fig. 6 is the amplitude, which is small. We'll see in the next section how this trans- 
lates into an expected AT in a given experiment, but clearly it will be difficult to detect 
these secondary anisotropics in CDM. Other cosmologies, in particular baryon isocurva- 
ture models, are more likely to produce a large secondary signal due to the Vishniac effect 
(Bond & Efstathiou 1987; Efstathiou 1988; Hu et al. 1993). Although we won't calculate 
the signal in these models here, we stress that our conclusion that Vishniac's term is the 
only important second order one applies in general. 
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7. Results and Discussion 



Now that we have all the physics under our belts, it is time to probe different ionization 
histories. In this section we present the C^'s for a variety of ionization histories and convolve 
them with the filter functions shown in Figure 1 to obtain the predicted signal in these 
experiments. 

FIG. 7. The visibility functions for the ionization histories we will discuss in this section. Shown for 
comparison are the standard recombination history and the fully ionized history. The small numbers 
alongside each line indicate the redshift of complete ionization. 

Figure 7 shows the visibility functions for the histories we will discuss. Note that the 
ionization in these particular histories is gradual so that even the history wherein ionization 
takes place at z = 10 differs from the standard recombination. The advantage of these 
particular histories is that they were generated in a consistent manner; i.e. a source of 
photons was injected continuously into the medium and the effects of recombination and 
electron heating were taken into account (Dodelson & Jubas 1992). The ionization history 
of the Universe could have been as sketched in any of the lines in Figure 7. The disadvantage 
of these histories is that in all of them re-ionization is gradual. We would like to be able 
to say something general about re-ionization without reference to these specific histories. 
Therefore, referring to them by the epoch of complete re-ionization is not a good idea. 

It is more useful to discuss the cumulative visibility function: J^" dr' g{TQ^ r'); this is 
shown in Figure 8. The cumulative visibility function is the probability that a photon has 
last scattered after a given time. Thus, for a no-rccombination Universe, Figure 8 shows 
that almost all photons scattered after r = .OStq. By contrast, the cumulative visibility 
function is only equal to .005 at r = .OStq for the standard recombination history. The 
other ionization histories lie somewhere in between. A nice feature of this number [the 
cumulative visibility function at r = .05ro] is that it characterizes each ionization history 
in an easily understandable way. Further, it lies between [standard recombination] and 
1 [no recombination] for all reasonable ionization histories. 

FIG. 8. The cumulative visibility functions for the ionization histories we will discuss in this section. The 
numbers by each denote the value of the cumulative visibility function when T = O.OStq. 

FIG. 9. The C/'s for different ionization histories ranging from standard recombination to no recombina- 
tion. Again the numbers give the value of the cumulative visibility function at T = .OStq. 
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For each of these histories we have calculated the expected signal in the CMB in the 
form of the C^'s. Figure 9 shows the results of the numerical integration of the Boltzmann 
equations. Working our way down from standard recombination, we see that if the cumu- 
lative visibility function is small at r = .05to [e.g. the curves labeled .10, .28, and .59], the 
effect of reionization is to damp out the primary anisotropics generated early on. The more 
effective the reionization [as parametrized by a larger cumulative visibility function], the 
more significant is the damping of the peak at I = 200. Indeed, if most photons scattered 
late [r > .OStq], then this peak goes away completely [e.g. the curves labeled .88 and .98], 
as we expect from our discussion in section 4 of the no recombination case. At the same 
time, though, secondary anisotropics - those generated at late times when velocities were 
larger - become important as reionization becomes more efficient. So for / < 50, Ci is 
actually larger in a Universe which reionizes early. 

We now convolve these C/'s with the filter functions for the experiments plotted in 
Figure 1. We are interested in the expected signal, as defined in Eq. (1.1). So we plot 
(ATg^^p^.)^/^ as a function of the reionization parameter. 

The expected signal in the Tenerife experiment is virtually independent of the ion- 
ization history. We see a small increase in the signal for no-recombination vs. standard 
recombination since the Tenerife filter samples part of the Doppler no-recombination peak. 
But the difference is very small, of order 5%, certainly too small to be meaningful at present. 

FIG. 10. The expected < AT^ >l/2 jj^ variety of experiments as a function of the reionization 
parameter, the cumulative visibility function at T = .05to. Small values of the cumulative visibility 
function correspond to standard recombination; values close to one correspond to no recombination. The 
signal in OVRO includes the second order Vishniac effect. 

The South Pole 91 filter is situated between the standard recombination peak at 
I ~ 200 and the no-recombination peak at / ~ 50. Since the former peak has a larger 
amplitude, the signal in SP91 is larger if the Universe had a standard recombination 
history than if it never recombined. It is interesting to note though that the curve is not 
monotonic: If the cumulative visibility function at r = .OStq lay between 0.5 and 1, the 
Doppler peak at I ~ 200 is significantly depleted and the secondary peak at / ~ 50 has 
not built up to its maximum value. Thus the expected signal in medium scale anisotropy 
experiments is lowest if roughly half the photons last scattered at the standard z ^ 1100 
and the other half at z ~ 100. The signal can drop by as much as a factor of two from 
what would be expected in standard recombination. Amusingly, the observed signal in all 
four channels was roughly a factor of two below what would be expected, but this must 
be taken with a grain of salt since there was probably a great deal of contamination from 
foreground sources (Gaier et al. 1992; Dodelson & Jubas 1993). 

The MAX experiment has a filter centered closer to the primary peak at / ~ 200. 
Therefore, the signal drops even more precipitously when the Universe is re-ionized and 
the primary peak is washed out. The signal can be a factor of three smaller than in standard 
recombination. Note again the effect of the secondary peak: as the cumulative visibility 
function reaches one, the secondary peak and, therefore, the signal increases. MAX and 
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SP91 are both very good probes of the ionization history, with the expected signal varying 
by a factor 0/ 2 — 3 depending on the ionization history. 

By far the best probe of ionization history is the Owens Valley Radio Observatory 
receiver. Note from Fig. 1 that the OVRO filter does not pick up the secondary peak at 
Z ~ 50. So the signal due to first order effects [pre-Vishniac] is completely negligible if the 
Universe never recombined. In fact, the expected signal drops by a factor of order ten if 
the Universe never recombined. This drop is large but it would be even more dramatic if 
not for the Vishniac effect, which produces a small, but non- negligible signal in the case 
of no-recombination when the primary signal has been completely washed out. 

To sum up, Fig. 10 gives meat to the common wisdom that the smallest scale exper- 
iments are most sensitive to the ionization history of the Universe. 

We started this investigation wondering whether a signal due to re-ionization could 
be misinterpreted as a primordial signal. We now know how a signal due to re-ionization 
differs from one due to standard CDM with standard recombination history, and we are 
confident that current and future experiments will probe such differences. What we have 
not done in this paper is explore how a signal due to re-ionization differs from a signal 
in other variants of CDM. Most troublesome are two variants most likely to be confused 
with re-ionization: (i) models where the primordial spectrum is not Harrison- Zel'dovich 
(n < 1) and (ii) models in which there are primordial tensor perturbations due to gravity 
waves (Crittenden et al. 1993). Both of these have the feature that the signal on large 
scales is the same as in standard CDM, while the signal on small scales is smaller than the 
standard one. They share this feature with re-ionized models. Therefore, distinguishing 
tilted or gravity-wave models from re-ionized models is a challenging task for cosmologists. 
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